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We generalize an ortho-normality relation between decay eigenmodes of equilibrium systems to 
nonequilibrium markovian generators which commute with their time-reversal. Viewing such modes 
as tangent vectors to the manifold of statistical ensembles, we relate the result to the choice of a 
coordinate patch which makes the Fisher-Rao metric euclidean at the invariant state, realizing a 
sort of statistical equivalence principle. Finally, we classify nonequilibrium systems according to 
their spectrum, arguing that a degenerate Fisher matrix is a signature of the insurgence of a class of 
phase transitions between nonequilibrium regimes. We exhibit an order parameter with power-law 
critical decay and prove divergent correlations between suitable unbiased estimators. 
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It is a chief objective of NonEquilibrium Statistical 
Mechanics (NESM) to draw lines and analogies between 
equilibrium and nonequilibrium systems, in search of as 
powerful instruments of analysis as are, for example, 
equilibrium ensembles. What makes a statistical ensem- 
ble equilibrium is certainly not its exponential form 

po<exp-£/3 Q X Q , (1) 

a 

nor the fact that it is a constrained maximizer of the 
Gibbs-Shannon entropy. Jaynes claimed that "essentially 
all of the known results of Statistical Mechanics, equilib- 
rium and nonequilibrium, are derivable consequences of 
this principle" pQ . Whether or not this program has been 
satisfactorily pursued, the maximum entropy inferential 
method might as well work for nonequilibrium steady 
states. Rather, nonequilibrium is discriminated by the 
fact that p is the terminal state of spontaneous relax- 
ation to a static scenario where no internal currents flow. 
The ensemble per se only tells a part of the story. The 
study of relaxation modes is essential to the full charac- 
terization of nonequilibrium behaviour. 

In this letter we estabilish a connection between ther- 
modynamical aspects of nonequilibrium systems, en- 
coded in their decay modes, and geometrical properties 
of the Fisher-Rao measure of distance between probabil- 
ity distributions. The former will be modelled through 
markovian evolution of a probability distribution p l on a 
finite sample space [2J, which decays exponentially fast 
towards an invariant nonequilibrium measure p. For such 
systems relative entropy 

Sip^^plHpl/P*)- (2) 

i 

plays the role of a Lyapunov functional. As was already 
known to Kullback [3] and thoroughly discussed by Baez 
[3] , relative entropy is also at the core of the definition of 
Fisher's matrix g a @ , which was reinterpreted by Rao as a 
metric g on the manifold of statistical states [5] . Already 
a standard tool in information theory and statistics [5J, 



the Fisher-Rao metric was rediscovered in equilibrium 
statistical mechanics [7], and more recently it is drawing 
attention in NESM iB] and in quantum information the- 
ory [9]. Singularities of <? Q/3 have been shown to signal 
the insurgence of quantum phase transitions [10] . 

Analogously, we will show that (i) the choice of a 
nonequilibrium generator picks out one preferred basis 
for the Fisher-Rao metric near the invariant state, in such 
a way that (ii) a diagonal Fisher matrix occurs for sys- 
tems which commute with their time reversal, which we 
refer to as p-normal, while (iii) degeneracies are accom- 
panied with the sort of critical behaviour which is typ- 
ical of nonequilibrium phase transitions — but with no 
cogent need for a thermodynamic limit on the number of 
states. The results are based on (iv) a generalization of 
the well-known [2j orthonormality relation between decay 
eigenmodes of detailed-balanced systems 

g a0 ._ y = §a0 (3) 

i Pi 

to systems whose generators commute with their time 
reversal. To this class there belong equilibrium systems, 
with real decay spectrum, and a class of p-normal, truly 
nonequilibrium systems, with oscillatory behaviour. We 
finally linger on the geometric meaning of the results, 
showing that (v) decay modes can be seen as local refer- 
ence frames for the Fisher-Rao metric, so that p-normal 
systems correspond to coordinate patches which make it 
euclidean at one point of the manifold, thus realizing a 
sort of statistical equivalence principle. As is customary 
in differential geometry, we will use a swarm of indices 
which are implicitly summed over when repeated. 

The system's finite sample space consists of n + 1 sam- 
ples, labelled by i,j,k. On it, a normalized probability 
density p l (a state of the system) undergoes markovian 
evolution, J^p* = Lp l . Conservation of probability re- 
quires the generator's columns to add to zero, 

L ' lJ ~ { - HkWki, i = j ^ 
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The u>'s are real positive transition rates, with units of 
an inverse time. Under mild assumptions on the con- 
nectedness of the state space and on the reversibility of 
paths, the markovian dynamics has a unique invariant 
state p, i.e. a null eigenvector of the generator, Lp = 0, 
towards which any initial state p e tends at large times, 
lim^oo e tL p e = p. The propagator exptL is a stochas- 
tic matrix with strictly positive entries. By a standard 
argument in the theory of Markov chains, the Perron- 
Frobenius theorem can be applied to prove that the n 
nonnull eigenvalues of L have negative real part, 



Lq a = (-l/T a + iu a )q a , 



l,...,n, 



(5) 



with T a a positive characteristic decay time and ui a a 
frequency. For the time being, we suppose the spectrum 
to be nondegenerate. The decay eigenmodes q a , plus 
the invariant state p, form a basis of eigenvectors of L. 
Propagating p e = EqP + e a q a up to time t yields 

-t/r a +iuj a t 



cxp(tL)p e = e p+ Y,' 



e a q 



(6) 



At late times decay modes are exponentially damped, 
so that by normalization e = 1, and consequently the 
eigenmodes' entries are seen to add to zero, qf = 0. 

With a slight abuse of language, we will say that L 
is of equilibrium when its invariant state satisfies de- 
tailed balance, with vanishing internal steady currents 
jij = WijPj - WjiPi, there being no net exchange of in- 
formation between samples. Following Andrieux [TT], we 
define the diagonal matrix ^Jp = diag{^/p 1 , . . . , -Jp A- 
Detailed balance holds if and only if matrix 



H:=-y/p- 1 Ly/p. 



(7) 



is symmetrical. Transformation ^ is a similarity of ma- 
trices. By the spectral theorem, H has a real spectrum 
which coincides with that of -L, He a = T ( ^ 1 e a , with ui a = 
and a ranging through 0, . . . , n. Eigenvector ef = ^/pl cor- 
responds to the steady state, with tq = oo. All others are 
related to the decay modes via ef = qfl^/pi. The spectral 
theorem also permits to choose orthonormal eigenvectors, 



uab. (a b\ cab 

h ■■= (e , e ) = o 



(8) 



where (•,•) denotes the euclidean scalar product. The 
l.h.s. provides the definition of the symmetric matrix 
h ab . Entries along the zeroth row yield the identities 
h oa = Y,iPi = 1 an d h 0a = Y,iQi = 0, which are true by 
construction also of nonequilibrium generators. Imposing 
h a/3 = 6 a P produces the desired result, Eq.|3]). Since, if 
we restrict to the field of reals, the derivation follows a 
chain of necessary and sufficient arguments, writing 



1 = { o g°e ) 



(9) 



we can state that g a ^ is diagonal in a basis of eigenmodes 
of a markovian generator with real non-degenerate spec- 
trum if and only if the latter is of equilibrium. 



On the complex field, the spectral theorem general- 
izes to normal matrices: it suffices (and is necessary) 
that HH^ = H^H to make H unitarily diagonalizable, 
i.e. with orthonormal eigenvectors with respect to the 
hermitian scalar product (v,w)c = Y,iV*Wi. Let us de- 
fine the time-reversal (also known as p-dual) markovian 
generator L = ^Jp 2 L 7 \/p 2 which also has p as its in- 
variant state [12]. Of course, there exists no generator 
which could possibly run the dissipative dynamics back- 
ward; L is what comes closest to it, as it inverts all of 
the steady currents ]Vj = -jji. Normality of H translates 
in the commutation relation of p-normal generators: 



[L,L]=0. 



(10) 



Since the entries of L are real, its 2k < n complex eigen- 
vectors come in complex conjugate pairs (<z+)* = qH, 
k = 1, . . . , k < n/2. It is then simple to see that Eq.( 10 ) 



is a sufficient and necessary condition for the real and 



imaginary parts of e£ = yjp 



9n 



ql to satisfy 



93 



(11) 



to be orthogonal among themselves, (SHe K , 3e K ) = 0, and 
to the remaining n-2k real eigenvectors e' , which satisfy 
Eq.|3]) on their own, (e l ,e' ) = 5 lL . 

Consider now the late-time behaviour of relative en- 
tropy with respect to the invariant state. We distinguish 
[n/2j phases in the space of generators with nondegener- 
ate spectra, parametrized by transition rates. 

(A) There are n real negative eigenvalues. Expanding 
Eq.([2]) to second order in e* near the steady state, with 
p 1 = ^/p(e° + e^e"), we obtain 



S-(p*b) 



I e a e 6\ e : e ) 



-\\ e* II 2 

2 II e Hg> 



(12) 



where first order contributions vanish. To second order, 
relative entropy is one-half the lenght of vector e t with 
respect to the positive-definite metric g al3 defined in the 
l.h.s. of Eq.([3|. Considering the explicit time evolution, 
Eq.|6]), with uj a = 0, we identify = e^ t l Ta e a . Nonequi- 
librium models with real spectrum then display superpo- 
sition of modes with different decay times, affecting the 
late time behaviour of relative entropy. This superpo- 
sition disappears for equilibrium generators. In a way, 
g a P measures the correlation between decay eigenmodes, 
with equilibrium modes being uncorrelated. 

(Bfc) There are n - 2k real eigenvalues -l/r t and k 
couples of complex conjugate eigenvalues —1/t k ± iu) K . 
In a basis of eigenmodes, the initial distribution reads 
p E = ^/p(e° + e L e L + e^e 1 ^), where a = ±. Notice that since 
p° is real one necessarily has (e~)* = e^. While relative 
entropy does not look like a real bilinear form when ex- 
pressed in terms of the complex vector e = (e^,e~,£ t ), 
it is indeed a positive bilinear form of the complexified 
vector £ = (D\e^ 7 3e^.,e L ). A direct calculation yields 

S(p*\ P) 



(13) 
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where the block-diagonal matrix 



n -\ e-^ ) 



(14) 



has k copies of the 2x2 rotation matrix R(tp), and the 
norm is to be calculated using the "complexified" metric 



\ 



r ) 



(15) 



Eq.(13l displays mixing of decay and oscillatory times, 
which only disappears when the generator commutes 
with its reversal and g a ^ is diagonal — with k blocks 
^ 1-2x2 ■ One last case is left out from our analysis: 

(C) The generator is defective, that is, degenerate 
eigenvalues lack a complete set of eigenvectors. This case 
and the phenomenology so far analyzed are better illus- 
trated with the aid of an example. Consider the following 
generator, parametrized by positive rates £, x 




X \ 
1 

■i-xj 



(16) 



The dynamics generated by L(£,x) is that of a hopping 
particle with a systematic bias in the counterclockwise 
direction, and one perturbed clockwise rate (see FigjT^,). 
Its phase space is depicted in Figfllj. By Kolmogorov's 
criterion, transition rates satisfy detailed balance if and 
only if the only macroscopic affinity vanishes, ln£/x 3 = 0, 
which traces the equilibrium line £ eq . The model corre- 
sponding to x = 1 = £ is known as the unbiased hopping 
particle, with twice degenerate eigenvalue A = 2 which 
affords a complete basis of eigenvectors. The space of 
parameters is partitioned into two phases of type A (in 
grey in Fig[TJi) and Bl (in white), marked out by the crit- 
ical lines £\ : \ + 3£ = 4 and £ 2 '■ £ = X- For the first class 
of models, direct calculation of the eigenvectors shows 
that g a P is diagonal only along the equilibrium line. In 
phase Bl one needs to turn to the complex components 
of the eigenmodes to be able to expand relative entropy 
as a positive bilinear form. Along I* are the biased hop- 
ping particle models L(l,x), which make the complexi- 
fied matrix diagonal. Their reversal if found by inverting 
the bias in the clockwise direction, yielding L = L T . 

With the exception of L (1,1), a generator L picked 
along the critical lines only has one eigenvector q relative 
to the degenerate eigenvalue -r _1 . A generalized eigen- 
vector u shall then be introduced, with Lu = -t^ 1 u + q, 
carrying L into Jordan's normal form. The time evolved 
exp(tL)p £ is seen to acquire a term <x te^ T q [TTj . 
Consider now a path T = {£(s), x( s )} m parameter space, 
as depicted in Fig[]J>. We first traverse the complex 
phase. At s = 1 we come upon an abrupt switch in the 
appropriate basis. Approaching the critical line l\ from 
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FIG. 1: a) Transition rates for the hopping biased particle 
with one clockwise perturbed rate, b) The parameter space 
of L(£,x): phases, critical and trivial lines, a path V. 

below, the imaginary part 3q+(s) becomes smaller and 
smaller until it vanishes; from above, the modes q+(s) 
and q~(s), respectively with higher and lower eigenvalue, 
tend to align. At the critical line the correlation matrix 
becomes degenerate and has a discontinuity 
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(17) 



An order parameter / s (oo) — of little physical meaning 
though — can also be engineered. Consider two vectors 
qt (s) and u 1 , respectively in the orthogonal complement 
of q~(s) and u. Projecting cxp[tL(s)]p £ -p along the 
two, and then taking the ratio, yields 

/ 2 > s >i(i) « [l+ce-*^- 1 -^ 1 )] -1 , f x (t) oc l/(l +c 'i). (18) 

The exponential decay, reaching an arbitrary nonzero 
value of f s (°°), becomes a power-law at the critical line, 
with / s (oo) = 0. Thence we took the liberty to refer to 
a class of nonequilibrium phase transitions, with an ac- 
ceptation that is reminescent of that employed for simple 
driven lattice models [T3]. This is also motivated by the 
equilibrium usage of the Fisher matrix, which coincides 
with the covariance matrix of the observable constraints 
X a which define the ensemble ([!]) . Crooks [S] commented 
that, as we vary the intensive parameters j3 a , correlations 
vary smoothly except at phase transitions, where diver- 
gencies occur. On this line, further insight might come 
from estimation theory. Consider an unbiased estimator 
e % = (e\, . . . , e l n ), whose average (e) Pe = e yields precisely 
the vector of parameters. The Cramer-Rao inequality 
[6], estabilishes a lower bound on the covariance matrix 
((£ - e)(e - e)) Pb > g^ 1 , where A> B means that A - B 
is positive semidefinite. Multiplying by g and taking the 
trace we obtain 



(lie- 



I sW 



(=) 



> n. 



s± 1,2. 



(19) 



A degenerate metric admits nonnull vectors of null norm. 
If these can be reproduced in terms of unbiased estima- 
tors, approaching the critical line degeneracy must be 
compansated by divergent correlations. 
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Indeed, a rich nonequilibrium phenomenology is 
marked by the peculiar representation of relative entropy 
near the invariant state, which can be interpreted as a 
metric on the manifold of statistical states V. Let us 
hint at its construction. Relative entropy is not a good 
distance: it is not symmetrical, and the triangle inequal- 
ity can be violated 0]. The way out of this puzzle is to 
stick to nearby distributions, as we did in Eq.(12), thus 



obtaining a local metric which measures the lenght of vec- 
tors e a q a living on the tangent space to V at p. When 
moving to a different neighbourhood, one will shift the 
reference probability distribution to p' , and there define 
the metric in terms of S(- \p'). If this procedure is carried 
on pointwise, one endows V with the Fisher-Rao metric. 
One can then assign coordinates x a to neighbourhoods 
of the manifold; associated to such coordinates is a basis 
of preferred tangent vectors d /dx a , which yield a matrix 
representative for the metric at each point of the neigh- 
bourhood. Notice that the Fisher-Rao metric is smoothly 
defined all over the manifold (except at boundaries and 
corners); it is its coordinatization that might suffer from 
pathologies, as is the case for our critical systems. 

But for n = 1, it can be shown that g has a nonnull Rie- 
mann curvature: while one can always choose a coordi- 
nate patch which trivializes the metric at one given point, 
there is no such coordinate transformation which simul- 
taneously makes g diagonal all over a neighbourhood. 
Given the twice-contravariant transformation law for the 
metric, g a '^' = A^'A^ g Q/3 , where A^' (i') = dx «/dx' a , is 
the inverse jacobian of the coordinate transformation 
x -*■ x (x), after Egs-ffl one realizes that the com- 
ponents of ef can be interpreted as the jacobian of an 
embedding patch, also called a frame, which trivializes 
the metric at p. Orthonormal frames are associated to 
p-normal systems; two such systems with the same in- 
variant state yield different orthonormal frames, which 
are connected by an 0(n + 1) "gauge" transformation. 

This mechanism lies at the heart of the Equivalence 
Principle of General Relativity. While gravity curves 
spacetime so as to prevent the definition of broad notions 
of "parallelism" and "simultaneity" , one can always find 
coordinates which make spacetime minkowskian at one 
point, and gravity indiscernible from a fictitious force. 
To a special observer, the frames' entries provide an in- 
ertial frame of coordinate axis: in a very precise way they 
measure how much the orientation of these axis differs, 
up to Lorentz transformations, from "bent" coordinate 
axis. This discrepancy is the gravitational field [H] . 

To conclude, the present letter estabilishes a close con- 
nection between the Fisher-Rao metric and markovian 
generators: The choice of a generator induces a natural 
identification of a positive semi-definite matrix represen- 
tation, nested in the late-time behaviour of relative en- 
tropy. A trivialized metric emerges for generators which 
commute with their time-reversal, including equilibrium 
systems; nonequilibrium phase transitions are induced 



by degenerate coordinatizations. This study calls for a 
careful treatment of the algebraic varieties of critical and 
trivial loci in the space of generators, and for a thermo- 
dynamical characterization of the order parameters. 

Normal systems seem to enjoy special properties, and 
deserve more in-depth study, in particular in relation to 
the fluctuations and large deviations of their microscopic 
jump trajectories. We notice at a passing glance that 
the "curvature" operator [L,L]5t 2 retains a geometrical 
flavour, as it measures how circuitation along infinitesi- 
mal parallelograms fails to close, when we first evolve the 
initial state with L up to time St and then run it back 
with L, rather then first evolving it back to -St with L 
and then run it forward with L. 

Finally, the neat framework suggests to deepen the 
nonequilibrium characterization of geometric objects 
such as Christoffcl coefficients, geodesic curves, intrinsic 
and estrinsic curvature. To conclude with a motto, we 
might claim that equilibrium systems are to nonequilib- 
rium thermodynamics what inertial frames are to gravity. 
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SUPPLEMENTARY MATERIAL 
On the Fisher-Rao metric 

The choice of a markovian generator L identifies a 
point on the manifold of statistical states and a set of 
decay modes. We first map the probability simpex into 
the surface of a sphere, then interpret the pushed- forward 
vectors' entries e? as the jacobian of a second coordinate 
transformation, in such a way that an equilibrium genera- 
tor corresponds to the choice of a coordinate patch which 
trivializes the Fisher-Rao metric at the invariant state. 
Given an invariant state, there exists a whole 0(n+ 1)- 
orbit of frame fields which trivialize the metric. Vicev- 
ersa, two equilibrium generators with the same invariant 
state yield gauge- equivalent frames. 

From a geometrical perspective — see for example [5]— 
the two matrices H and —L are both representations of 
an operator H with eigenvectors 



„a d 



'* dzi ■ 



(20) 



The peculiar notation employed for the basis vectors de- 
notes that yjp might be seen as the inverse jacobian of 
the coordinate transformation pj h» Zj (p) = which 
maps the probability simplex {pi e [0, l] n+1 : T,jPj = 1} 
into a portion of the hypersphere with square radius 
Y,i zf = 4. Both are embeddings of the abstract mani- 
fold V of probability distributions into R™ . Each such 
coordinatization (say, x) of "positions" on V endows the 
(n + 1) -dimensional vector space V = E Tl+1 of "veloci- 
ties" , attached to p, with a preferred basis of directions 
a /dxi. Vectors e Q span the n-dimensional tangent space 
T p V c V, while e° describes how a neighbourhood of p 
sits in the embedding space. 

A metric at p is a positive semidefinite bilinear form 
h p [v, w], with vectors v, w e V attached to p. A choice of 
coordinates endows the dual space V * with a basis of lin- 
ear forms dx a , such that dx a [ 9 /dx b ] = S b . Vectors' compo- 
nents are then obtained by projecting v a = dx [v], while 
the metric in coordinates reads h p = h ab (x)dx a <X> dxb- 
Again, the notation "dx a " highlights the transformation 



properties of linear forms: if x <-* x'(x) is a diffeomor- 
phism, the metric's entries transform twice contravari- 
antly (dropping dependencies) 



,a 6 



dxg dx b 
dx' a , dx' b , ' 



h ab 



(21) 



A metric is fully characterized by its action on a complete 
set of vectors. In our case, we define h p via 



h p [e a ,e b ] := ft <%e) 



--: S tj e, <g> e 



[ dx a ' dx b ] ( 22 ) 



where Eq.([8]) is inbuilt in the l.h.s., while the r.h.s. helps 
introducing the frame forms e, = efdx a BacJ. Given a 
metric, frames are determined up to 0(n + l)-rotations 
>-* R^Bj, which leave the metric's coefficients unal- 
tered. When such rotations are performed point-by-point 
all over a neighbourhood of the manifold, in a differen- 
tiable way, one talks of a gauge transformation. 

Comparison of Eq.(22) with Eq.(21) reveals that ef 



is the inverse jacobian of a coordinate tranformation 
z x(z) which pulls back the metric's components to the 
unit matrix. For example, choosing "spherical" coordi- 
nates x a = 8 l a Zi, one obtains the euclidean lenght element 
h v = (eiz)® 2 . The intrinsic metric at p, g p , is defined as 
the restriction of h p over the hypersurface V, acting on 
the tangent space T V V . By construction e° is orthogo- 
nal to all other e Q 's. Finally, the Fisher- Rao metric on 
the manifold is simply the disjoint collection of such local 
intrinsic metrics, g = \_\ P nV gp- 



On relative entropy in the complex phase 



Consider the probability density 

p E = yfp [e° + e L e L + § (e£ + ie 2 K ) (ej - ie K 2 ) + c.c] , (23) 

where el = SKe£, e£ = 3e+ = 3e*. Here k ranges among 
the k couples of complex conjugate eigenvectors, I spans 
the remaining n-2k real modes. Propagating up to time 
t yields, in full extent, 



-t/r L 



e L e L + Y, e t/T " [ cosiuj^elel - sm(u; K t)E* e£ + sin(w K t)e 2 e K x + cos(oj K t)e 2 K e^] 



(24) 



In short, e tL p, 



*yp~(e° + e T Q}e), where e = (e^e^e 1 ), 
i = (el ,££,£') and Q* is defined in Eq.([l4|. Formula (fl3]) 



then easily follows. Notice that when L is p-normal, from 
the hermitian orthogonality relations 

(ei + ie 2 ,ei + ie 2 )c = 0, (ei - ie 2 , e\ + ie 2 )c = c (25) 



there follow the real orthonormality relations 

(ei,e 2 ) = 0, (ei,ei) = (e 2 ,e 2 ) = f. (26) 

Whatever the normalization chosen for the complex 
eigenvectors, real and imaginary parts have the same nor- 
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malization, so that g consists of blocks of type cl2x2- It 
follows that for p-normal systems oscillatory late-time be- 



haviour in the relative entropy, encoded in ft^gft 1 , will 
disappear, since R{w K t) T ■ cl 2x 2 ■ R(u K t) = cl 2x 2- 



On the three-state system 



The off-diagonal element of the complexified Fisher ma- 
trix reads 



9 12 



(x-1) 



(x 4 + 2x 2 y + Ax 3 y + x 4 y + 6x 2 y 2 + 7x 3 y 2 + 



P1P2P3 

3xy 3 + I3x 2 y 3 + x 3 y 3 + 6xy 4 + 5x 2 y 4 - y 5 + 6xy 5 ) . 



We work out in full extent the three-state model, with 
x = £iU = X- Equilibrium holds when the Kolmogorov 
criterion is satisfied, that is then the ratio of products of 
clockwise over counterclockwise rates yields 1. Its locus 
defines the equilibrium line £ eq : x - y 3 = 0. The charac- 
teristic polynomial of L is 

det(Al-L) = A 3 + (x + 2y + 2>)\ 2 + (3y 2 + 2x + 2y + xy+l)\, 



with roots A = and A ± = - (rr. + 3y + 2 ± \/A) /2 where 
the discriminant A = (x - y)(x + 3y - 4) vanishes at the 
critical lines i\ : x + 3y - 4 = and £2 '■ x - y = 0. The 
unnormalized invariant state is 

p=(y 2 + y+l,y 2 + y + x,y 2 + xy + x). (27) 

We proceed analyzing the three cases of interest. 

(A: A > 0). The eigenvectors are 

' 4y - (x + y ± \/A) 2 \ 
q ± = -Ax + 2y(x + y± \/A) 
k 2x 2 + 2xy - 4y 2 ± 2x\/A J 

The off-diagonal element of the Fisher matrix reads 

+- = [4y - (x + y + 7A) 2 ][4y - (x + y - x/A) 2 ] 

y 2 + y + 1 

[-4a; + 2y(x + y + VA)] [-Ax + 2y(x + y- \/A)] 
y 2 + y + x 

[2x 2 + 2xy - Ay 2 + 2x\fK] [2x 2 + 2xy - Ay 2 - 2x\fK] 
y 2 + xy + x 

Imposing g 12 = 0, we find the solutions: x - y 3 = 0; 
x 2 - y 3 = 0; 1 + 2x + 2y + xy + 3y 2 = 0. The first is precisely 
the equilibrium line, which corroborates our conclusion. 
The third equation has no physical solution. The sec- 
ond solution is more tricky, as we didn't expect its ap- 
pearence. There the metric is not diagonal, as q ± are not 
the eigenvectors of L(x, x 2 ! 3 ). In fact in this case q+ = 0, 
and there exists a correct complete basis of eigenvectors: 



/ 


2-2/3 


) 




1 


-x 2 ! 3 - 


X \ 




-x 2 / 3 - 


X 


q- = 




x 2 ' 3 




V 


X 


j 




V 


X 


1 



It vanishes on the line £* : x = 1, but there is another 
class of solutions in the first quadrant. It can be shown 
that those solutions correspond, as above, to systems for 
which q ± arc not the correct eigenvectors. The time- 
reversal operator is 



L(x,y) 



-x-y 

x+y+y 2 
1+y+y 2 

xy+xy +y 
1+y+y 2 



y+y +y 

x+y+y 2 

-1-y 

x+xy+y 2 
x+y+y 2 



x+xy+xy 2 y 
x+xy+y 2 

%y+y 2 +y 3 

x+xy+y 2 

-i-y / 



(28) 



(29) 



Setting the commutator with [L, L] to zero, yields 



[L,L]= {X - I){y3 - X) M = 

P1-P2- P3 



with 

Mu = -y(x + 2xy + 2y 2 + y 3 ) 

M12 = (l + y + y 2 )(x + xy + 2y 2 ) 

M13 = -(x + y 2 )(l + y + y 2 ) 



We don't need to specify more entries of M, as system 
(291 has as its unique solutions x = y 3 (the equilibrium 
line) and £* : x = 1, as expected. 

(C: A = 0). Notice that on the critical lines q + = g_ 
and 3q + = 0. Looking for eigenvectors, on £2 we obtain 



Al - L(x, x) 




(30) 



whose kernel is 1-dimensional, unless x = 1. Similar con- 
clusions can be drawn on i\. 



(B: A < 0) Eigenvectors q± have real and imaginary 
parts given by 

' 2x - x 2 + y 2 - 2xy ^ 
D\q + = y 2 + xy - 2x , Jq + 



K x + xy - 2y 



-x-y 

y 

x 
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